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Abstract 

The aim of this paper is to prove existence of weak solutions of hyperbolic-parabolic 
evolution inclusions defined on Lipschitz domains with mixed boundary conditions de¬ 
scribing, for instance, damage processes and elasticity with inertia terms. To this end, a 
suitable weak formulation to deal with such evolution inclusions in a non-smooth setting 
is presented. Then, existence of weak solutions is proven by utilizing time-discretization, 
iJ^-regularization of the displacement variable and variational techniques from [HKll] to 
recover the subgradients after the limit passages. 

AMS Subject classifications 35L20, 35L51, 35K85, 35K55, 49J40, 49S05, 74A45, 74G25, 
34A12, 35K92, 35K35; 

Keywords: Hyperbolic-parabolic systems, doubly nonlinear differential inclusions, existence 
results, energetic solutions, weak solutions, linear elasticity, rate-dependent damage systems. 


1 Introduction 

The gradient-of-damage model motivated by Premond and Nedjar in [FN96] describes the 
damage progression by microscopic motions in solid structures resulting from the growth of 
microcracks and microvoids. In this approach, an internal variable z models the degree of 
damage in every material point. It is bounded in the unit interval [0,1] with the following 
interpretation: the value 1 stands for no damage, a value between 0 and 1 qualifies partial 
damage and the value 0 indicates maximal damage. Beyond that, elastic deformations are 
described by a vector-valued function u which specifies the displacement from a prescribed ref¬ 
erence configuration D. The evolution law for u and consists of two equations: a hyperbolic 
equation for the mechanical forces and a parabolic equation for the damage process involving 
two subgradients. The time evolution of {u, z) under constant temperature can be deduced 
from the laws of thermodynamics and is summarized in the following hyperbolic-parabolic 
PDE system (see [FN96, Frel2]): 

utt - div (We(e(u), z)) = £, (la) 
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Zt - ApZ + VF, 2 (e(n), z) + f'{z) + ^ + 99 = 0, 
where the subgradients ^ and are given by 


i e 5 /[o,oo)(2) 




If G 5/(_oo,0](2:t) 


with the snbdifferential 

{ 0 } if 2 ; > 0 , 

(— 00 , 0 ] if 2 ; = 0 , 

J ifz< 0 , 

with the subdifferential 

{ 0 } if Zt < 0 , 

[ 0 ,oo) if 2 :f = 0 , 

0 if -St > 0 . 


oo,0] (^t) 


The system is snpplemented with the following initial-boundary conditions: 

u = b on Td X (0, T), 

W^e{^{u), z) ■ V = ^ on Tn X (0, T), 

V 2 : ■ v = 0 on do,, 

tt(0) = in n, 

ut(0) = in n, 

z(0) = z° in n. 


(lb) 


(2a) 


(2b) 


(3a) 

(3b) 

(3c) 

(3d) 

(3e) 

(3f) 


The hyperbolic eqnation (la) is the balance eqnation of forces containing inertial effects 
modeled by uu, the parabolic eqnation (lb) describes the evolution law for the damage 
processes and (2a) as well as (2b) are snbgradients corresponding to the constraints that the 
damage is non-negative {z > 0) and irreversible {zt < 0). The symbol ApZ := div(|V|^“^V 2 :) 
denotes the p-Laplacian of 2 ;. The inclnsions (2a)-(2b) are explained in more detail below. 
Moreover, i denotes the exterior volume forces, / a given damage-dependent potential, e{u) 
describes the linearized strain tensor, i.e. e(u) = ^(Vu -|- (Vu)"*"), and Td and Tn indicate 
the Dirichlet part and the Neumann part of the boundary 912. The elastic energy density W 
is assumed to be of the form 

W{e, z) = ^h{z)Ce : e, (4) 

where C is the stiffness tensor and h models the influence of the damage. We assume h' > 0 
and that complete damage does not occur, i.e., h is bounded from below by a positive constant. 
Let us give an interpretation for the system (lb), (2a) and (2b) modeling damage processes 
via gradient flows (for a physical motivation by means of microscopic force balance laws and 
constitutive relations we refer to [FN96]): For this purpose, we observe that the inclusion 
(2b) is equivalent to the complementarity problem 


Zt < 0 , (fzt = 0 , ip>0. 


( 5 ) 


Now, we introduce the free energy T" as 


J^iu,z) := + W{e{u),z) + /(^) + L[o,oo)(^)^ dx, 


( 6 ) 
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where the indicator function /[o,oo) can be interpreted as an obstacle potential, i.e., the damage 
variable z is forced to be non-negative. The gradient-of-damage term models influence 

of the damage on its surrounding and also has a regularizing effect from the mathematical 
point of view. A calculation reveals </? = —zt — C with C, G dzF{u, z), where tp is given by (2b) 
and C is a subgradient of the generalized subdifferential dzJ-{u,z), i.e. 

C =+ W^zi£iu),z) + f'{z) + G 9 -^[ 0 ,oo)(^)- (7) 

The complementarity formulation implies 

if Zt < 0 then zt G —dzT{u, z), 

< if Zt = 0 then C < 0 ; 

Zt > 0 not allowed. 

We gain the following interpretation of a non-smooth evolution: As long as the driving force 
—C given in (7) is non-positive, the evolution is described by the gradient flow zt = —C tending 
to minimize the free energy. Whenever ^ becomes negative, zt is 0. We remark that also an 
activation threshold for the damage process can be incorporated by adding a linear term to 
the potential /. 

The aim of this paper is to give a notion of solution to the system (l)-(3) considered on 
bounded Lipschitz domains 0 in a weak sense and to prove existence of weak solutions. In 
the following, we summarize the main difficulties we were faced with and how they are solved: 

• Let us point out that even on smooth domains severe difficulties arise in establishing 
strong solutions to (l)-(3). Because due to a missing viscosity term —div(V(z)e(ut)) in 
(la), no L^(77^)-regularity estimates for ut are available. This, in turn, leads to severe 
difficulties in obtaining higher regularity estimates for the damage variable z. In fact, to 
obtain L°°{L?‘)-a. priori estimates for the p-Laplacian ApZ and for the subgradients (2a)- 
(2b), we may consider Moreau-Yosida type regularizations {Cs^'^s) for (^, 99 ). Then, we 
may test (lb) with {—ApZ+^siz))t, integrate and use the estimates {ips{zt), —ApZt) > 0 
and {ifsizt), i^siz))t) > 0 and (—ApZ,^ 5 (z)) > 0. But to conclude the desired estimates 
via integration by parts, we would require L'^-bounds for the strain rate s{ut) (cf. [RR14, 
Seventh a priori estimate]). 

In our case the situation is even worse since no global 77^-elliptic regularity results are 
available for bounded Lipschitz domains and mixed boundary conditions. 

In consequence, we have to devise a concept of weak solutions. In the paper [HKll] 
(where the inertia term uu is neglected), a notion was introduced, which allows us to 
formulate the double inclusion in system (lb), ( 2 a) and ( 2 b) with fairly weaker regu¬ 
larity assumptions. The notion combines a variational approach with a total energy- 
dissipation inequality. Here, we adapt this formulation to the present situation since 
the inertia term uu in (la) leads to new terms in the energy-dissipation inequality. 

• To prove existence of weak solutions in the sense mentioned above, we are required to es¬ 
tablish a total energy-dissipation inequality. However, when using a time-discretization 
scheme, the discretized system may exhibit error terms which converge to 0 in the 
time-continuous limit if L^-bounds for the strain tensor e(u) are available. 

To this end, we hrstly study a regularized version of system (l)-(3) where an additional 
fourth-order term gives naturally rise to spatial L7^-regularity for the displacement 
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variable rt in a weak setting. Finally, we perform a limit analysis to obtain a weak 
solution of the system (l)-(3). 

• In order to perform the limit passage of the subgradients occurring in the two approx¬ 
imations (namely the time-discretization and the Ff^-regularization), approximation 
techniques from [HKll] (see Lemma 2.1 and Lemma 2.2) are utilized. Although the 
precise structure of the subgradient ^ can be deduced in the Ff^-regularized version 
of (l)-(3) (see (28)), the situation is more involved when the passage from the 
regularized to the limit system is investigated. In contrast to [HKll] and by the best 
knowledge of the authors, no strong L^(FI^)-convergence of u can be established in this 
limit passage due to the additional inertia term. 

We solved this problem by exploiting the precise structure of S, in the regularized system 
such that it cancels out with other terms in (lb) on certain parts of the domain (see 
(70)). For the remaining terms we apply lower-semicontinuity arguments and uniform 
convergence of the damage variable to pass to the limit system. Only then we are able 
to recover the subgradient 

We would like to conclude this introduction by comparing our contribution with existing 
mathematical works involving gradient-of-damage models and some of their regularization 
strategies in the literature (the choice of literature is, of course, only an excerpt and not 
complete): 

• Among the pioneering works, [FKS99] (see also [FKNS98] for quasi-static evolutions) 
analyzed the damage model in [FN96] in the case of an one-dimensional rod. The 
authors considered an inertia term in the force balance equation and a viscosity term 
in the stress tensor. They were able to prove local-in-time existence of solutions. In 
the case of reversible damage processes, uniqueness of solutions was proven. Several 
different approximation strategies came into play such as penalization methods for 
the subgradients, a regularization of the damage rate function, truncation and time- 
redarding methods. 

• A 3D elasticity-damage model with inertial effects was analytically studied with a reg¬ 
ularization of the damage rate function in [BS04] and with a viscosity term in the stress 
tensor in [BSS05]. Local-in-time existence was proven by means of Schauder fixed-point 
theorem and uniqueness was obtained in the case that either irreversibility or bound¬ 
edness of the damage variable is dropped. Homogeneous Dirichlet boundary conditions 
for the displacements and smooth domains were assumed. 

• Rate-independent damage models with a quasi-static balance of forces were investigated 
in [MR06]. This work covers existence results for the rigid body impact problem. By 
means of a so-called energetic formulation and the usage of F-convergence, the case of 
complete material desintegration (complete damage) is also discussed. 

• The work [MTIO] analyzed further rate-independent damage models and proved ex¬ 
istence of energetic solutions by a new method for the construction of suitable joint 
recovery sequences. Also temporal regularity properties of the solutions were shown. 

• The transition between rate-dependent and rate-independent damage models, the so- 
called vanishing viscosity limit, was investigated in [KRZ13a]. The authors assumed 
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a quasi-static equilibrium of forces and a higher order Laplacian in the damage law 
(for 3D). Existence of energetic solutions was shown, where the proof is based on a 
regularization of the damage rate function. Furthermore, a priori estimates were derived 
and uniqueness results were deduced for special cases. 

• A coupled system of damage, viscoelasticity with inertia and heat conduction was stud¬ 
ied in [RR14]. In case of isothermal processes, smooth domains and homogeneous 
Dirichlet boundary conditions for the displacements, existence results were established 
and uniqueness was proven for reversible damage processes. 

• [HKll] as well as [HK13] explored PDF systems coupling damage processes with quasi¬ 
static elastic systems and Cahn-Hilliard equations for phase separation. A notion of 
weak solutions involving variational inequalities and a total energy-dissipation inequal¬ 
ity was introduced. The authors prove existence of weak solutions and higher integra- 
bility of the strain tensors. They also provide abstract approximation techniques (see 
Lemma 2.1 and Lemma 2.2) which will be utilized in the present work to study (l)-(3) 
involving an inertia term in the force balance equation. In particular, this system allows 
for the presence of elastic waves interacting with the damage evolution. 

Let us also mention that the p-Laplacian —ApZ := —div(|Vz|^“^Vz) in (lb) describes 
the diffusion of damaged material parts across their proximities. For mathematical 
reasons p will be chosen to be larger than the space dimension because the embedding 
^ C'(D) plays an essential role in the mentioned approximation lemmas. By 
replacing the operator —Ap by —{6Ap + A) in the damage law and performing a limit 
analysis <5 0, we were even able to handle damage models with the standard Laplacian 

(i.e. p = 2 as done in [HK13]). However, to keep this presentation short, we will not 
follow this approach here. 

Structure of the paper 

In Section 2, we introduce some notation and preliminary mathematical results from [HKll, 
HK13]. The main part is Section 3. We state and justify a notion of weak solutions in 
Subsection 3.1. The proof of the existence theorem ranges from Subsection 3.2 to Subsection 
3.3. At first, we prove existence of weak solutions for an iL^-regularized problem by using 
a time-discretization scheme and by applying variational techniques from Section 2 to pass 
to the time-continuous system. Finally, we get rid of the regularization by a further limit 
passage which is performed in Subsection 3.3. 

2 Notation and preliminaries 

Throughout this work, let p € (n, oo) be a constant and p' = pl{p — 1) its dual and let 
D C (n = 1,2,3) be a bounded Lipschitz domain. For the Dirichlet boundary Td and 
the Neumann boundary Tn of 5D, we adopt the assumptions from [Berll], i.e., Td and Tn 
are non-empty and relatively open sets in (9D with finitely many path-connected components 
such that Td n Tn = 0 and Td U Tn = d^l. 

The considered time interval is denoted by [0, T] and := D x [0, t] for t G [0, T]. Form 
now on the time-derivative in front of the main variables and the data is denoted by dt and 
the partial derivative of the density function W (see (4)) with respect to one of its variables 
e and z is abbreviated by W^e and respectively. Furthermore, we define for k > 1 the 
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spaces 


:= {u e W’^’P{Q) I u > 0 a.e. in Q}, 

W’:’P{n) := {m e W^’P{n) I M < 0 a.e. in Q}, 

:= |m S H^{Q) I m = 0 on Fd in the sense of traces}. 

The following variational and approximation results are crucial to establish existence of weak 
solutions. They are aimed to pass to the limit in certain integral inequalities with weakly 
convergent functions constraining the set of test-functions (see “Step 1” and “Step 2” in the 
proof of Proposition 3.8). Roughly speaking, Lemma 2.1 provides an approximation result 
in a strong topology for test-functions satisfying weakly convergent constraints. This enables 
us to pass to the limit in the corresponding integral inequalities. Then, Lemma 2.2 provides 
a method to drop the limit constraints on the set of test-functions, where, in turn, a new 
integral term arises on the right-hand side. This method will be used to perform the limit 
passage of a time-discretized version of (lb) and to recover the subgradients (2a) and (2b). 
We will frequently make use of the compact embedding 

w^’P(L!) cm 

without mentioning in our considerations. In the following, the notation {C = 0} 3 {/ = 0} 
for functions in L°°{0,T;W^’P{Q)) should be read as 

{x G n I C(x, t) = O} 3 {x G n I /(x, t) = O} for a.e. t G (0, T). (8) 

Beyond that, the subscript r always refers to a sequence with \ 0 as fe ^ oo. 

Lemma 2.1 (See [HK13]) Let 

• frj eL°°{0,T-,wl’P{Q)),T>0 

with frit) —)• f{t) weakly in IT^’^(n) as t \0 for a.e. t G (0,T), 

• C e L°°{0,T;Wl’P{Q)) with {( = 0} D {f = 0}. 

Then, there exist a sequence fr £ L°°(0, T; iy_}’^(r2)) and constants > 0 such that 

• Ct ^ C strongly in L‘i{Q,T-,W^^P{Ll)) as r \ 0 for all q G [l,oo), 

• Cr —>■ C weakly-star in L°°{0,T;W^’P{LI)) as t \ 0 , 

• Cr < C o,.e. in LIt for all t > 0 (in particular {Ct = 0} 3 {C = 0} ), 

• W,tCr(i) < frit) in LI for a.e. t G (0, T) and for all r > 0. 

If, in addition, C ^ f n.e. in LIt then the last condition can be refined to 

Ct < /r n.e. in LIt for all r > 0. 

Lemma 2.2 (See [HK13]) Let f G LP'{Ll;W^), g G L\Ll) and z G W];P{Ll) with f-Vz > 0 
a.e. in and {/ = 0} 5 {z = 0} in an a.e. sense. Furthermore, we assume that 

[ (/ • VC + dC) dx > 0 for all f G W^’P{Ll) with {C = 0} 3 {z = 0}. 

Jn 

Then 

[ (/ • VC -|- gC) dx > f max{0, (7}C dx for all C G W^’P{Ll). 

Jq J { 2 = 0 } 
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Remark 2.3 In [HK13], g is assumed to be in L^[Q). But the proof extends to g £ 
without any modifications. 

3 Analysis of the hyperbolic-parabolic system 

3.1 Notion of weak solutions and existence theorem 
Assumptions on the coefficients 

Besides the general setting introduced at the beginning of the previous section, we now state 
the conditions needed for the coefficient functions in the PDE system (la)-(lb) in order to 
develop a weak notion and to prove existence results: 

We assume / G ^^([0,1],M+) for the damage potential function (occurring in (lb)) and W 
to be given by (4) with h G ^^([0,1];M) and h > g on [0, 1] for a constant g > 0. Moreover, 
we will use the assumption h' > 0 on [0,1]. The stiffness tensor C should satisfy the usual 
symmetry and coercivity assumptions, i.e. Cijik = Cjuk = Cikij and e : Ce > colep for all 
e G and constant cq > 0. 

Weak formulation 

The main idea for a weak formulation is to rewrite the doubly nonlinear differential inclusion 
in system (lb)-(2b) into a variational inequality and a total energy inequality. This kind of 
notion was introduced in [HKll] and is adapted to the present situation. 

Lemma 3.1 Let the data 

G G L^(fl;M"'), G with 0 < < 1 a.e. in 11, 

£ G T2(o,r;L2(51;M’^)), 
b G H\0, T; M”)) n ^^(0, T; M”)) 

and the functions 

u G L2(0, T; K")) n r1(0, T; R^(0; M")) n ^^(O, T; L^{LL- M")), 

z G TP(0, T; n 1T^’P(0, T; IT^’P(fl)), 

^,y,eL\0,T-,L\n)) 

be given. The following statements are equivalent: 

(i) The functions {u,z,^,ip) satisfy the PDE system (l)-(3) pointwise in an a.e. sense to 
the data ,v^, z^ ,i,b). 

(a) The functions {u, z,(^,(p) satisfy the following properties: 

• initial and Dirichlet boundary conditions: 

u{0) = u^, dtu{0) = v^, 2:(0) = z^, u = b onTj) X (0,T), (9) 

• force balance: 

{dttu.OHi + f We(e(u),2) : e(C)dx = f £-Cdx (10) 

Jn Jn 

for all C G and for a.e. t € (0,T), 
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damage evolution law: 


0 — y (^{dtz)( + |V 2 ;|P • VC + W^zi^iu), z)C + f'iz)C^ dx + (Cj Oiyi.p + Ov^i.i 

( 11 ) 

for all C £ and for a.e. t G (0,T), 

conditions for the subgradients and the damage variable: 

for all C G a.e. t£ (0,T), (12a) 

a.e. in Q, X {0,T) (12b) 


dtz < 0 


and 


for all C G a.e. t £ (0,T), 

a.e. in Q X (0,T), 


z{t))wrp < 0 
z>0 

total energy-dissipation balance: 

T{t) + lC{t) + V{0, t) = T{0) + /C(0) + Wext{d, t) 
for a.e. t G (0,T) with 

■= y + bb {e{u{t)),z{t)) + fiz{t))^ dx, 


(13a) 

(13b) 


(14) 


^(i) := J ^\dtu{t)fdx, 

rt 


free energy: 
kinetic energy: 

dissipated energy: P(0, i):= / / \dtz\^dxds, 

Jo Jn 

WextiO,t) := [ [ W^e{£{u),z) : e{dtb)dxds - [ [ dtu - dubdxds 

Jo Jn Jo Jn 

+ f dtu{t) ■ dtb{t) dx — f dtu{0) ■ dtb{0) dx 

J Q J n 


external work: 


+ 


f f i ■ dt{u — b)dxds. 

Jo Jn 


Remark 3.2 The essential feature in the notion presented in (ii) is that we only require 
((/j,C)vKi’P < 0 in (12) instead of the full variational inequality ((/?, C — dtz)\Yi,p < 0. In 
view of the complementarity formulation (5) it means that the condition {q:>,dtz)yyi,p = 0 is 
dropped. However, as we will see in the proof, the total energy-dissipation balance allows us 
to recover this identity. 


Remark 3.3 Let us give a physical interpretation of the term bVea;t(0,t). To this end, re¬ 
member that the Cauchy stress tensor is given by 

a = W^e{s{u),z). 
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Now, by using integration by parts, the force balance equation (la) as well as the boundary 
conditions in (3), Wexti^N) transforms into 

yVext{0,t) = f f {a ■ dtb) ■ udxds + f [ ( - diY{a) ■ dtb + duu ■ dtb) dx ds 
Jo Jtu Jo Jn ^ ' 

+ f [ I ■ dt{u — h) dx ds 

Jo Jn 

n [a ■ u) ■ dtb dx ds + / / £-dtudxds. 

The first integral term may be interpreted as the work performed by the prescribed (time- 
dependent) Dirichlet boundary data and the second integral term is the work performed by the 
external forces £. 

Proof of Lemma 3.1 

To “(ii) implies (i)”: 

The equations (la) and (lb) can be recovered by standard arguments from (10) and (11). 
The inclusion (2a) follows by the variational inequality (13). 

It remains to prove the validity of the inclusion (2b). As mentioned in Remark 3.2, we need 
to show {ip,dtz)iYi,p = 0. 

To this end, denote the free energy functional without the indicator part (see (6)) by 

T{u,z) := J ^\Vz\^-\-W{e{u),z)-\-f{z)^ dx. 

The Gateaux derivatives duJ^ and d^J-” are given as follows: 

{du^iu,z),0Hi = f W^ei£iu),z) : £{C)dx, 

Jn 

{d,T{u,z)X)w^,p = [ {\Vz\P-^Vz-VC + W,Mu),z)C + f'{z)C) dx. 

Jn 

Testing (10) with dtu — dtb (note that dtu — dtb = 0 on Td) yields 

{duT{u,z),dtu)Hi = {du^{u,z),dtb)Hi + f i-dt{u-b)dx- f ^hdtul'^ dx 

Jn Jn z 


+ {dttu,dtb)Hi 

Testing (11) with dtz shows 

{dzF{u,z),dtz)wi,p = - 


Jn 


dx - {f,,dtz)wi,p - {ip,dtz)iYi,p. 


(15) 


(16) 


We obtain by adding (15) and (16), integrating over time and using the chain rule for 
ft^iu{t),z{t)): 

F{u{t), z{t)) - T{u{0), 2 :( 0 )) = J (^{duT{u{s), z{s)), dtu{s))Hi + {d^T{u{s), z{s)), dtz{s))wi,p'j ds 

= - / ]-\dtu(t)\'^ dx + [ ^|5iu(0)|^dx + f [ W^e{£:{u),z) : £{dtb)dxds 
Jn z Jn z Jq Jq 


9 



+ 


[ [ i ■ dt{u — b)dxds+f f dttu ■ dtbdx ds 

JQ Jn Jo Jn 

{Cidtz)w^’P ~ 



Furthermore, integration by parts in time shows 

nt c pt 


f f dttu ■ dtb dx ds = — f f dtu ■ dub dx ds + [ dtu{t) ■ dtb{t) dx — f dtu{0) ■ dtb{0) dx. 
Jo Jn Jo Jn Jn Jn 

(18) 


By applying (18) to (17), we obtain 

J~{t) — J~(0) = —IC{t) + /C(0) + yVext{d, t) — 11(0, t) — dtz)\yl,p — {(f, dtz)\yl,p. 
Subtracting the total energy-dissipation balance (14) from above yields 

{Cj9tz)y^i,p + {<f,dtz)iYi,p = 0. (19) 


In order to obtain {ip,dtz)y/i,p = 0, we will prove (^,(9t2;)^i,p = 0 and then applying (19). 
Since ^ G L^(0, T; L^(n)) by assumption, we find by (13) 

/ ?(C “ - 2 ;) dx < 0 

Jn 

for all C G F^(n) and for a.e. t G (0,T). This shows ^ = 0 a.e. in {z > 0} and ^ < 0 a.e. in 
{z = 0}. Since dtz = 0 holds a.e. in {z = 0}, we obtain 

{i,dtz)wi,p = / idtzdx = d. 

Jn 

To “(i) implies (ii)”: 

The force balance equation (10), the damage evolution law (11) and the variational inequality 

(13) follow from (la),(lb), (2a) and the boundary conditions (3a)-(3c) without much effort. 
It remains to show the variational property (12) and the total energy-dissipation balance 

(14) . 

In fact, (2b) implies the complementarity formulation 


dtZ < 0, {^P,dtz)y/l,p — 0, (<y9, > 0 (20) 

for all C £ W^^{VL). In particular, (12) is shown. 

To show (14), we consider the calculation (17) which follows with the same arguments as 
in (i). Taking into account the conditions {ip,dtz)-^Y\,p = 0 (see (20)) and {^,dtz)y^ri,p = 0 
(follows as in (i)), we obtain the total energy-dissipation balance. □ 

From Lemma 3.1 we make the following observations: 

• The weak formulation in Lemma 3.1 (ii) requires much less regularity as assumed there. 
Since we avoid the full variational inequality {ip, C — dtz)-[Yi,p < 0 by means of the total 
energy-dissipation balance, the notion makes still sense for z G L°°(0, T; kF^’^(n)) H 
77^(0, T; L^(n)). This is very important because the a priori estimates will give no 
better regularity in this weak setting (remember that is only a bounded Lipschitz 
domain). 


^ dxds 
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The regularity u e L°°(0, T; (0; T; W^))nH^{0, T; M”))*) 

for u is sufficient for the notion in Lemma 3.1 (ii). We observe that the term VVext avoids 
the second time-derivative of u. In fact, we may also write (by using (18)) 


>Vext(0,t) = f [ We(e(tt),z) : e((9t6) dxds + f {dttu{s), dtb{s))ds 
Jo Jn Jo 

+ [ [ i ■ dt{u — b) dx ds 

Jo Jn 


provided that u € {Q]W^))*) (note that dtb is not necessarily 0 on Td). 

• Since ^ is assumed to be in L^(0,T; L^(n)), the W^’^-dual product rewrites as 


= f C{t)Cdx. 

Jn 

The regularity assumptions for ^ can be weakened to L^{0,T; (Q)) since the integral 

term on the right-hand side exists in this case for all C € VL^’^(n) provided that p £ 
(re, oo) by employing the embedding W^’P{Q) ^ L°°{Q). 

• Conditions (11) and (12a) can be reformulated as one inequality as 

0 < [idtz)C + iV^r^Vz • VC + LL..(e(n), z)C + /'(^)c) dx + (C, C)w^,p (21) 

holding for all C S W^^{Vl) and for a.e. t G (0, T). This eliminates cp in the weak 
formulation. 


A further not so obvious observation is stated in the following lemma: 

Lemma 3.4 Let u, z, C cmd p as well as u^, , z^, I and b be as in Lemma 3.1. Suppose 

that the eonditions (ii) in Lemma 3.1 except the total energy-dissipation balance are satisfied. 
Instead we assume for a.e. t £ (0,T) the total energy-dissipation inequality 

T{t) + Kit) + P(0, t) < Fid) + /C(0) + Wextid, t), (22) 


where the funetions F{-), /C(0, ■), T>{d, •) and VV’ea;t(0, •) are defined in Lemma 3.1 (ii). 
Then, the total energy-dissipation balanee is satisfied, i.e., (22) is an equality. 


Proof. We find (17) and {f,,dtz)\yi,p = 0 with the same argumentation as in the proof of 
Lemma 3.1 (see “(ii) implies (i)”). 

Then, taking 

{ip,dtz)iYi,p < 0 


into account (which follows from (12) tested with p = dtz), we obtain the desired “>”-part 
of the total energy-dissipation balance. □ 


These observations motivate the following notion of weak solutions which is equivalent to the 
PDE system (l)-(3) provided sufficient regularity. 
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Definition 3.5 (Notion of weak solutions) Let the data {u^,v^, ,i,b) be given as in 

Lemma 3.1. A weak solution of the PDE system (l)-(3) is a triple of functions 

u G L°°(0, T; {Q; M")) n T; L‘^{Q; R^)) n H^{0, T; (O; M"))*), 

2 G L°°{0, T; W^'P{Ll)) n H\0, T; L‘^{LI)), 

^ G L^{n X (o,r)) 

satisfying (9), (10), (12b), (13), (21) and the total energy-dissipation inequality 

F{t) + lC{t) + P(0, t) < T{0) + /C(0) + Wext{0, t) (23) 


for a.e. t £ (0, T). 

The main aim of this work is to prove existence of weak solutions in the sense of Definition 

3.5. 

Theorem 3.6 (Existence of weak solutions) Let the assumptions in Section 2 and on 
the coefficients f, h and C be satisfied. Furthermore, let the data ,v^, z^ ,£,b) according 
to Lemma 3.1 be given. Then, there exists a weak solution of system (l)-(3) in the sense of 
Definition 3.5. 

We emphasize that proving the total energy-dissipation balance (14) in the weak setting 
seems to be out of reach by the authors’ best knowledge. However, by utilizing non-local or 
regularized versions of the Laplacian operator in the damage law (lb), the energy balance 
can be recovered as shown in [KRZ13a, KRZlSb]. 

3.2 Existence of weak solutions for an iT^-regularized system 

We firstly study a regularized version of the PDE system (l)-(3) in order to prove Theorem 

3.6. The enhanced regularity allows us to control an error term which occurs in the discrete 
version of the energy-dissipation inequality. The passage to the limit system is then performed 
in the next subsection. 

The regularized PDE system is described in a classical notion by a quadruple {u, z, ip) of 
functions satisfying the following equations pointwise: 


dttu— (i\v{W^e{^{u),z)) -b 5div( div(V(Vu))) = 1, 

(24a) 

dtz - ApZ -b W^zi£iu),z) + f'{z) -b ? + 93 

= 0 , 

(24b) 

f, G ^-^[ 0 , 00 ) (^)) 


(24c) 

p G dfi_oofi]{dtz) 


(24d) 

with the initial-boundary conditions 

u = b 

on Pd X ( 0 ,r), 

(25a) 

(We(e(u),z) -5div(V(Vu))) • = 0 

on Pn X (0,T), 

(25b) 

(5V(Vtt) • = 0 

on dLl X (0, T), 

(25c) 

0 

II 

> 

on dLl, 

(25d) 

u( 0 ) = dtufif) = 2 ;( 0 ) = z^ 

in Ll. 

(25e) 
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In order to treat the fourth order regularization term “(5div( div(V(Vu)))” with the given 
constant 5 > 0 analytically, we introduce the linear operator A : —>■ 

by 

{Au,v)h2 := / (V(Vu), V(Vu))Knxnxn dx := ^ / dxiXjUkdxiXjVkdx. 

According to the regularization we modify the weak formulation in Definition 3.5 as follows: 

Definition 3.7 (Notion of weak solutions for the regularized system) We consider the 
given data as in Lemma 3.1. Additionally, we assume A 

weak solution of the regularized PDE system (24)-(25) is a triple of functions 

u G L°°(0, T; M”)) n IDb°°(0, T; L^{VL- M")) n T; (D; R^))*), 
z G L°°(0, T; iybP(!^)) n T; L^{n)), 

^ G L\n X (o,r)) 

satisfying (9), (12b), (13), (21), the regularized version of the forces balance equation 

{dttuX)m + [ W^e{£{u),z) : e{C)dx + 6 {Au,C)h^ = [ ^'Cdx (26) 

Jn Jq 

for all ( G //p^(D;M"’) and a.e. t G (0,T) and the regularized version of the total energy- 
dissipation inequality 


Tsit) + lC{t) +V{Q,t) < + 1C{D) + Wixt{d,t) (27) 

for a.e. t G (0,T), where T, 1C, V and Wext CLre given as in Definition 3.1 and 

Fs{t) := F{t) + -{Au{t),u{t))H 2 , 

:=Wext{0,t)+6 f {Au{t),dtb{t))H2dt. 

Jo 

The goal of this section is to prove the following proposition: 

Proposition 3.8 (Existence of weak solutions for the regularized system) Let the as¬ 
sumptions in Section 2 and at the beginning of Section 3.1 on the coefficients f, h and 
C be satisfied. Furthermore, let 6 > 0 and the data {u^,v^, z^ ,£,b) according to Lem¬ 
ma 3.1 be given. Additionally, assume G H‘^{Q]R^), £ G (^^’^(O, T; L^(D; M”)) and 
b G C‘^d{0,T; H'^{Q-,R'^)). Then, there exists a weak solution of system (24)-(25) in the 
sense of Definition 3. 7. 

Moreover, for this weak solution, the subgradient ^ has the form 

f. = -X{2=o}max|o,lT,^(e(u),z) + /'(z)|. (28) 

We will complete the proof of Proposition 3.8 at the end of this section. The proof is based 
on a time-discretization scheme. To this end, let {0, r, 2r,... , T} be an equidistant partition 
with fineness t :=T/M > 0 of the interval [0, T]. The final time index is denoted by M G N. 
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By a recursive minimization procedure starting from the initial values := (u^,z^) 

and u~^ := — tv^, we obtain functions {u^^, z'^^) for every m = 0,..., M. To this end, we 

hx an m G {1,..., M} and define the functional M”) x —)• M by 


:= (^^\Vz\P + W{e{u),z) + fiz)-i'^-u^ dx + ^(^'u,'u) 




+ 


r 

^ ^T,S 

+ y 


2 

T 

r2 



L2 



L2 


A minimizer of in the subspace 

X Z -5 := [u G | u\r^ = 6™)|r„} x {z G W^’P{n) \0<z< z™-'} 

obtained by the direct method is denoted by The velocity field is dehned as 

the time-discrete derivative and the discretizations 6™ and are set to 


■.= b{mT), i^:=i{mT). 


(29) 


For a discretization re™ G 6™}, we introduce the piecewise constant inter¬ 

polations w, w~ and the linear interpolation w with respect to the time variable as 


w{t) := re™ with m = [t/r] , 

w-{t) := with m = (t/rl , 

w{t) := + (1 - ^ ^ ^ ^ ^ t-{m-l)T 

T 

and the piecewise constant functions tr and t~ as 

tr '■= \t/T~\ T = min{mr | m G No and mr > t}, 
t~ := max{0,tT- — r}. 


We would like to remark that by above definition we have w(t) = w{tr) for all t G [0, T]. 
Within this notation, the velocity field satishes 


dtVr{t) 




with m = [f/r] and t G (0,T], 

Since (5 > 0 is assumed to be constant in this section, we mostly omit the subscript <5 in 

m m j -pm pm 

We obtain the following time-discrete (in)equalities by the minimizing property of the func¬ 
tions (tt™,z™) with respect to the functional W™ over x ZJ^. 


Lemma 3.9 The functions («(", z™) G x Zft^ for m = 0,..., M and, consequently, the 
piecewise constant and linear interpolations 


Ur,Vr G L°°(0,T;/72(ff;]R")), Ur,% G W^’^{0,T-,H^{n;R^)), 
Zr G L^{0,T-W^’P{n)), Zr G W^’°°{0,T-W^’P{n)) 


satisfy 


14 














(30) 


(i) for all C G and for all t G (0,r); 

/ dtVr-Cdx+ / W^e{£{Ur),Zr) ■■ e{C)dx + 5{AUr,OH2 = / ir'C^X, 

Ju Jo. JO. 

(a) for all C G and all t G (0,T) with 0 < C < ®-e- in 

0 < ^ (jVZrr^VZr • V(C - Zr) + {W,MUr).Zr) + f'{Zr) + dtZr){f - Zr)) dx. (31) 


Proof. The minimizer (tt™,z™) fulfills the variational property 

-d„j-™«, 2 r) = o, 

where Nz^^zlf) denotes the normal cone to Z'f^ at z™. Equivalently, 
n™ - 2u^-^ + n!"-2 


/ 


2 

f'o ' 


Cdx+ [ lT,e(e(n™), 2 :™) : e(C) dx + 6{Aufd, C)m = / • C dx 

J Q, J Q, 


for all C £ 

0 < ^ • v(c - c(i)) + (iE,,(e(n™),2r)+ /'(^d + 


_ ^m—1 


)(C-^. 


for all C G 


□ 


Lemma 3.10 (A priori estimates) 

(i) The following a priori estimates hold uniformly in r and 6: 

V^lkT,<5||L°°(0,r;_ff2(O;R")) < ll_H'2 + C, 

lkT,<5||L°°(0,T;//l(f^;®"))nWl.°°(0,T;L2(0;R")) < '^C\\u^\\h2 +(7, 

lbT,5||L°o(0,T;I,2(0;Rn)) < y/5C\\u^\\H2 + C, 

lbT,<5|lL°o(0,T;L2(0;R"))nHl(0,r;(H2^(0;R'i))*) ^ '^C\\u^\\h2 +C. 


(32a) 

(32b) 

(32c) 

(32d) 

(32e) 


(ii) For fixed <5 > 0, the following additional a priori estimates hold uniformly in t: 

lkT,(5||L°°(o,T;tvi'P(r2)) ^ C*, (33a) 

lkT,5||L°°(0,T;tVl'P(r2))nfl'l(0,T;L2(O)) < C. (33b) 

Remark 3.11 The a priori estimates in (ii) are due to the -regularization for u (the a 
priori bound (32a) to be more precise). Later on, it will also enable us to establish the total 
energy-dissipation inequality in the limit regime r \ 0. Then, the total energy-dissipation 
inequality will give a priori estimates of type (ii) uniformly in 5. 


dx. 
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Proof of Lemma 3.10 

To (i): Proof of the a priori estimates (32a)-(32e): 

Testing (30) with Ur — u~ — {hr — h~) and using the estimate 


/ 


dtVr ■ {Ur — Ur)dx > 


|2 

Il2 


- - 


2 

L2 


as well as the convexity estimates (note that Zr > Zr) 

[ W^e{£{uT),ZT) ■ £{ur - Ur)dx > [ [W{e{ur)■, Zr) - W {e{u~), z~)) dx (34) 
Jn Jn 

(35) 


yield 

1 

2 

+ 


< 


s s 

S{AUr,Ur - U~)h2 > -{AUr,Ur)H2 - -{AUr , U~) h2 , 


1 (5 S 

L2 - 2 hr(.t)\\l2 + -{AUr{t),Ur{t))H2 - -{AUr{t),Ur{t))H2 

{W{s{Ur{t)),Zr{t)) -W{s{Ur{t)),Zr{t))) dx - [ dtVr{t) ■ {br{t) - b~ {t)) dx 

Jn 

[ irit) ■ {Ur{t) - Ur{t) - {br{t) - br (t))) dx 

Jn 

+ [ W^ei£{Urit)),Zr{t)) : £{br{t) - b~{t))dx + 6{AUrit),brit) - br{t))H 2 - (36) 

Jn 


The right-hand side can be estimated by Young’s inequality as follows (r/ > 0 denotes a freely 
chosen constant) 


r.h.S. < CT [\\ ir { t )\\ j^2 + \\ v ^ 


2 -L 
1,2 + 




^2 + 'n\\^{uT{t))\\l2 + Cr,\\£{dtbr{t))\\‘j^2 


-h 5{AUr{t),Ur{t)) h2 + 6{Adtbr{t),dtbr{t))H2 


(37) 


Summing (36) over the discrete time points r,2r ,... ,tr, using the estimate (37), an Lp‘{L?‘)-a 
priori bound for ir and an H^{H^)-a. priori bound for br, we obtain for small rj > 0 


1 „ 

- \\Vi 
2 " 


2 m 

L2 + 


^{AUr{t),Ur{t))H2 + c\\£{Ur{t))\\\2 - [ [ dtVr ■ dtbr dx ds 

^ Jo Jn 

^ \ lk°llL2 + + J^W{£{u^),Z^)dx 

+ C f {l + \\Vr{s)\\‘l2 +r]\\£{Ur{t))\\‘l2 + 6{AUr{s),Ur{s))H2 +Crj) ds ( 38 ) 

Jo 

for all t G [0,T]. The discrete integration by parts formula yields for all t G [0,T] 

f f dfVr ■ dfbr dx ds = f Vr{t) ■ dtbr{t) dx — [ V^-dthr{d)dx 

Jo Jn Jn Jn 

_ £’ r . dMs) - aM> - r) pg) 

Jo Jn T 
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Applying (31^ to (38) and using an a priori bound for the second discrete time- 

derivative of br, we eventually obtain for all t € [0, T] 


\ Il^r(i)|li2 ^{AUr{t),Urit))H2 "h c||e(u^ (t)) ||^2 

rU 


< -h < 5 ||u °||^2 + ^ (Ibr(s)|li2 + r]\\s{Ur{s))\\l2 + 6{AUr{s),Ur{s))H2) ds 
/■*"/■ ^ dtbr{s) - dtbr{s - t) 

/ / [s) ■ 

Jo Jq 


dx ds 


< -I-( 5 ||t (°||^2 -I- J (||ur(s)||^2-I-||u^ (s)||^2-I-?/||e(t(,-(s))||^2-I-( 5 (Aut-(s), Ut-(s))j^ 2) ds 

We conclude by a discrete version of Gronwall’s inequality (32a)-(32d) and 

l|W,5||L°°(0,T;L2(n;M")) < 'J'^C\\u^\\h2 +C 

with the help of Korn’s inequality. By using these a priori estimates, a comparison argument 


in (30) shows 


sup {dtVrit),C) 


dt 


< 


5||Ur(t)||^2(o,T;//2(n;R")) “I" f j (l ^,e(^('*^r) j ) 1 + |^t| ^ dx dt 


< (5C||u' 


0||2 

1^2 


+ c. 


Thus (32e) is proven. 

To (ii): Proof of the a priori estimates (33a) and (33b): 

Testing (31) with Zr which is possible due to 0 < Zr < a.e. in 12 x (0, T) and using the 
I • l^-convexity estimate 


yield 


IVZrP ^VZr • V(Zr — Z:_ ) dx > -IIVz. 

P 


P _l||Vz“ll^ 

T\\LP pW^’^rULP 


VZr{t)\\l^ - ^I|VZ^ {t)\fLP +r\\dt%{t)\\l 2 

< [ {W^zieiUrit)), Zr{t)) + f'{Zr{t))) {z~ {t) - Zr{t)) dx. 


(40) 


Thus, by Young’s inequality [r] > 0 denotes a freely chosen constant) 

l||V2,(()|li, - ^\\Vz;(t)\\l.+r\\a,Z{mh 

< Tr]\\dtZr{t)\\l 2 +TCrj\\W^z{e{Ur{t)),Zr{t)) + f'{Zr{t))\\l 2 . 


Summing over the discrete time points r, 2 t, ... ,tr and using the constraint 0 < Zt- < 1 a.e. 
in 12 X (0,T) as well as the a priori bound Il'Wr||L°°(o,T;VKi’4(n;R")) < C (which follows from 
(32a) for fixed <5 > 0 and the continuous embedding i7^(12;M’^) iy^’'’^(12;M”)), we end up 

with (33a) and (33b). □ 

The a priori estimates in Lemma 3.10 give rise to the following convergence properties which 
follows from standard weak and Aubin-Lions type compactness results [Sim86]. 
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Lemma 3.12 There exist functions 


u € L°°(0, T; M”)) n T; L'^{n-, M")) n H^{0, T; (Jl; M"))*), 

z e L°°(o, T; w^’P{n)) n h\o, t- L^{n)) 


satisfying (9), (12b) and (13b) and a subsequence \ 0 as k co such that 



weakly-star in L°°(0, T; M"')), 

(41a) 


strongly in L'^(0, T; iy^’^(n; M”)) for every q G [l,oo), 1 < s < 2*, 

(41b) 

Ur^ u 

weakly-star in L°° (0, T; n (0, T; (52; M*")), 

(41c) 


strongly in L'^(0,T; 1T^’^(52;M”)) for every q G [l,oo), 1 < s < 2*, 

(41d) 

'^Tk , v~ dtu 

weakly-star in L°°(0, T; T^(n; M”)), 

(41e) 

Vr^ dtU 

weakly-star in L°°(0, T; L^(n; M”)) and weakly in H^{0,T] (i/p^^ (fl; M”))* 

(41f) 


weakly-star in L°°(0,T; 1T^’^(52)), 

(41g) 


strongly in L'^(0,T;L°°(52)) for every q G [l,oo), 

(41h) 

Zrk Z 

weakly-star in L°°(0,T; 1T^’^(52)) and weakly in 2/^(0,T;L^(52)), 

(41i) 


strongly in C{Q x [0,T]) 

(41j) 


as k oo for fixed 6 > 0. The constant 2* denotes the critical Sobolev exponent. Moreover, 
we obtain the following convergence properties of the data 



strongly in L^(0,T;L^(52;M"")), 

(42a) 

Kk -t b 

strongly in 2/^(0, T; i7^(52; M”)), 

(42b) 

(Jtbrk ~ dtbrki' ~ J'k) ^ Q ^ 

strongly in L^(0, T; ff^(52; M"")). 

(42c) 


Tk 


Proof. To (42a)-(42c): 

We set X := H‘^{Q-,W^). By exploiting the fundamental theorem of calculus for functions 
with values in X and the assumed Lipschitz continuity of dub in time, a straightforward 
calculation shows 


f 


dtbrkit) - dtKkit - Tk) ^ 

2 „ 

f [ idttb{i) - dttb{f))dids 

Tk 

X Jo 

1 

1 

1 t- 


dt 


X 


-A / (f f ||9it6(0 - 9tt6(t)||^dids^ dt 

^k JO \JTf. Js-Tk J 


T / rt 


1 

^k Jo \Jt 




C — t| di ds ) dt 




= 4(7V|r ^ 0 


'S-Tk 


<2Tk 


as k CO. Thus (42' 


) ID OllUW ii. 
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To (41a)-(41j): 

Standard weak and weak-star compactness resnlts applied to the a priori estimates in Lemma 
3.10 reveal existence of functions 

u, u- G 

u G L'^(0, T- H^{VL- M”)) n T; 

v, v- G L°°(0,r;L2(L!;M”)), 

V G L°°(0,r;L2(L!;]R")) n (L!;M"))*), 

G L°°(0,r;lT^’P(L!)), 

T G L“(0, T- kT^’P(L!)) n F^(0, T; L^{a)) 

satisfying (9), (12b) and (13b) and subsequences indexed by such that 

Ur^. —>■ u weakly-star in L°°(0, T; Lf^(n; M”)), 

^ u~ weakly-star in L°°(0, T; Lf^(n; M”)), 

Ur, weakly-star in 

Vr, —>■ V weakly-star in L°°(0, T; L^(fl; M"')), 

v~ weakly-star in L°°(0, T; L^(n; M"')), 

Vr, —>■ n weakly-star in L°°(0, T; L^(fl; M"")) and weakly in Lf^(0,T; (fl; M”))*), 

Zr* ^ ^ weakly-star in L°°(0, T; IT^’^(n)), 

^ weakly-star in L°°(0, T; IT^’^(fl)), 

£ri; ^ ^ weakly-star in L°°(0,T; IT^’^(n)) and weakly in 

as k oo. Taking into account 

Ur, - Ur, = ndtUr, 0 strougly in L°°(0, T; L^(0;M"')), 

we obtain u = u~ = u. Analogously, we get v = v~ = v and z = z~ = z. The identity 
dtUr, = Vr, implies dtu = v. 

Therefore, we obtain 

u G L°°(0, T; M”)) n 1T^’°°(0, T; ^^(Ll; M”)) n H^{0, T; (L!; M”))*), 
z G L°°(0, T; lyi’P(Ll)) n H^{0, T; L2(L1)) 

such that (41a), (41c), (41e), (41f), (41g) and (41i) is satisfied for the subsequence {Tk}kefi- 
Compactness arguments (in particular using the embedding (7(0) valid for 

p G (n,oo) and Aubin-Lions type results; see [Sim86]) show (41b), (41d), (41h) and (41j). □ 


Remark 3.13 By choosing further subsequences (we omit the additional subscript), we also 
obtain for fixed 5 > 0 


Vur,,Vuf^,Vu 


u 




Vu 


^r,,Zr,,Zr, ^ ^ 


pointwise a.e. in Q x (0,T), 
pointwise a.e. infix (0,T), 
pointwise a.e. infix (0,T) 
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and for a.e. t G (0, T) 


weakly in 
weakly in 


u{t) 
^ z{t) 


as k oo. 

Strong convergence of V^r^. Vz in L'P{Q x (0,T)) can be shown by a subtle approximation 
argument introduced in [HKll], 

Lemma 3.14 There exists a subsequence of Tk (omitting the additional subscript) such that 
Zr^ in LP{0,T; ILi’P(fl)) ask ^oo. 

Proof. We apply the preliminary result cited in Lemma 2.1 based on the work [HKll]. By 
Remark 3.13, the sequence fulfills the assumption. According to Lemma 2.1 there 

exists an approximation sequence {Cr*.} ^ with the properties 

Crj, ^ 2 strongly in L^(0,T; (43) 

0 < Crfe < pointwise a.e. in H x (0, T) for all A: G N. (44) 


We omit the subscript k for notational convenience. Property (44) enables us to test (31) 
with ^T-. Integration over the time variable shows 

[ [ \VZr\^~‘^VZr-V{Zr-CT)dxdt < f ( (W^zi^iUr), Zr) + f (Zr) + dtZr) {(t - Zr) dx dt. 

Jo Jn Jo Jn 

A uniform p-monotonicity argument and the above estimate show (c > 0 is a constant) 


c||V 2 — Vz. 


LP(0,T;LP(C;R")) 


< 


10 Jn 

rT 


{\Vz\P-^Vz - \VZr\P~'^VZr) ' V(z - Zr) dx dt 


= [ [ \\/z\P-^Vz-V{z- Zr)dxdt+ f [ \VZrf-^VZr-V{Zr-CT)dxdt 

Jo Jn Jo Jn 

+ ! [ \^Zr\P-^VZr-V{Cr-z)dxdt 

Jo Jn 

< f [ {W^zi£iUr),Zr) + fiZr) + dt%) iCr - Zr)dxdt 

Jo Jn 

+ [ [ \Vz\P-^Vz-V{z-Zr)dxdt+ [ [ \Vzr\P-^Vzr-V{Cr-z)dxdt. (45) 

Jo Jn Jo Jn 

In the following, we prove that every term on the right hand side converges to 0 as r 0. 

• The first integral on the r.h.s of (45) can be estimated as follows 

[ [ {W^z{e{Ur),Zr) + f'iZr) + dtZr) (Ct - Zr)dxdt 

Jo Jn 

^ II ld^,2(^(Ur)) ^r) + / (^t) ||x,2(-o^T;L1(0)) 11^"^ ~ ||L2(q^2’;L°°(0)) 

+ ||'9i2T||^2(Qy)||CT - ^T||^2(f^j,)- (46) 
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By using the boundedness of Ur in L°°(0, T; M"’)), boundedness of dt'Zr in L^(0, T; L‘^{Q,)) 

(see Lemma 3.10), boundedness Zr G [0,1] a.e. in Qt and the convergence properties 
(41h) and (43), we obtain convergence to 0 as of the two summands on the right hand 
side of (46). 

• Due to the convergence (41g), the second integral on the r.h.s. of (45) converges to 0 
as T \ 0. 

• We estimate the third integral on the r.h.s. of (45) by Holder’s inequality: 

f j • V(Ct ~ 2:) dxdt < ~ 2 )||lp(0,T;Lp(S1))' 

Because of the boundedness property (33a) and the strong convergence property (43), 
we obtain convergence to 0 of the integral term above. 

Combing the convergence result Vzr Vz strongly in ^^(0, T; LP(D; M”)) as r \ 0 with 
Lemma 3.12, the claim follows. □ 


The notion of weak solutions as given in Definition 3.7 requires the validity of the total 
energy-dissipation inequality. However, in this discrete setting, we are only able to prove 
an approximate version of this inequality. But the 77^-regularization enables us to recover 
the postulated total energy-dissipation inequality in the limit r \ 0 as already indicated in 
Remark 3.11. 


Lemma 3.15 For a.e. t G (0,T) the approximate energy-dissipation inequality 

J^,{t) + JCrit) +VriO,t) + £r{0,t) < .T(O) +/C(0) + (0, t) 

with 

Fr{t) := j (^^\VZr{t)\^ + W{£{Ur{t)),Zrit)) + f{Zr{t))^ dx-\-^{AUr{t),Ur{t))H^, 


(47) 


^r{t) := f ^|up(t)pdx, 

Jn ^ 

Pt-( 0,t) := / / \dtZr\^dxds, 

Jo Jn 

rt 


VVea:t(0>*) := / [ W^e{£{Ur), Zr) : £{dtbr) dx ds - [ [ 

Jo Jn Jo Jn 


r 


-I- f Vr{t) ■ dtbrit) dx — f ■ dtbr{0) dx d- f [ £r ■ (dtUr — dtbr') dxds 

Jn Jn Jo Jn ^ ' 

-\-6 {Aur{s),dtbr{s))H 2 ds 

Jo 


and the “error term 

£. 


(o,t):= r [ 

Jo Jn 


1 h{z.r ) — h{Zr) 


Ce{u.^ ) : £{u.^ ) dx ds 


+ / / W^z{£{Ut) Ft) dtZr dx ds 

Jo Jn 
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ri 


f{Zr) - f{z^ ) 


dxds+ / f\zr) dtZr dx ds 

Jo Jq 


holds. 


Proof. By employing the estimate (which is slightly sharper than the convexity estimate (34)) 

W^e{^{Ur), Zr) '. e{Ur — Uf) dx 
> [ {W{e{ur),Zr) -W{s{uf),zf)) dx 

Jn 


I 

Jn 


1 


+ / o (^(^r ) - H^t)) Ce('u^ ) : e('u^ ) dx, 


we obtain by testing (30) with Ur — uf — {br — bf) (cf. (36)): 

1|I /mi2 1|I-/\I|2 d , , ,,, 5 

- ki. 

2 " 


L2 - 2 11^^ + 2^^Urit),Ur{t))H2 - (0)//2 

+ [ {W{s{Ur{t)),Zrit)) -W{s{uf{t)),Z~{t))) dx - [ dtVr{t) ■ {br{t) - bf {t)) dx 
J Q, J Q, 


^ Ini “ Hzrit))) C£{u^ (t)) : e(«^ (t)) dx 


< 


[ ir{t) ■ {Ur{t) - (t) - {br{t) - 6^ (t))) dx 

Jn 

+ f W^e{£{Ur{t)),Zrit)):£{br{t)-bf{t))dx + 6{AUr{t),brit)-bf{t))H2. (48) 
Jn 

By testing (31) with Zr, we obtain 

^l|V2:r(i)|liP - ^||Vz“(t)||^p +T||at£r(t)|li2 

< [ {W^zi£{Ur{t)),Zr{t)) + f'{Zr{t))) {Z~{t) - Zr{t))dx. (49) 

Jn 

as in the proof of Lemma 3.10. Adding the estimates (48) and (49), summing over the discrete 
time points and taking into account formula (39) yields (47). □ 

We are now in the position to establish the equalities and inequalities of the weak formulation 
of Dehnition 3.7 by passing r \ 0. As before, we omit the subscript k in the sequence {Tk}k&N- 

Proof of Proposition 3.8 

The functions u and 2 : from Lemma 3.12 already satisfy (9), (12b) and (13b). It remains to 
show (13a), (21), (26) and (27) from Definition 3.7. 

To (26): We find for all C G 


I 

Jn 


dtVr{t) • C dx = {dtVr{t), C)m 


by using the canonical embedding L^(D;M”) “—)• (L7p^(D;M"))*. Keeping this identity 
in mind, integrating (30) over time from t = 0 to t = T and passing to the limit 
r 0 for a subsequence by using the convergence properties in Lemma 3.12, we obtain 
a time-integrated vesion of (26). Then, switching back to an “a.e. in t”-formulation 
shows (26). 
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To (13a) and (21): The limit analysis for these equations are performed in two steps and 
makes use of the approximation technique cited in Lemma 2.1 and the extension result 
cited in Lemma 2.2. 

Step 1: Let C € with {C = 0} T {z = 0} (see (8)). By Lemma 2.1, we 

obtain a sequence {Cr^} ^ (we omit k) and constants > 0 

with the properties: 

Cr ^ C strongly in L^(0,T;lT^’^(fl)), (50a) 

0 > i^r,tCr(i) > —Zrit) in 0 for a.e. t S (0, T) and all r > 0. (50b) 

The property (50b) and Zr < z~ holding pointwise a.e. in 12 x (0, T), we also find 
for a.e. t € (0, T) 

0 < ^'T-,tCr(i) + Zr{t) < z~{t) a.e. in 12. 

In consequence, for a.e. t € (0, T), we can test (31) with + Zr{t) and obtain 

^T,t 

> 0 . 

We divide this inequality by the positive constant and integrate over the time 
interval [0,T]. The time-integration is necessary to exploit the weak convergence 
property for dt'Zrit)) in L^(0, T; L^(12)). More precisely, we use the convergence 
properties in Lemma 3.12, Remark (3.13) and Lemma 3.14 to pass to the limit 
T \ 0 for a subsequence and end up with 

0 < / [ {\Vz\P-‘^Vz ■ VC + (W^(e(u), z) + f'iz) + dtz)C) dx dt. 

Jo Jn 

In particular, we get an a.e. in time t formulation. 

Step 2: We may apply Lemma 2.2 to the above variational inequality. Then, we obtain 
for all C £ IT'1’^(12) the inequality 

"'-I Li dtzC + \Vzr^Vz-VC+{W,Mu),z) + fiz) + ^C) dxdt (51) 

with C £ L^(12 X (0,T)) given by 

-X{2=o}max|o,(9tz -k W^^{e{u),z) + f'{z)y 

Due to dtz < 0 a.e. in Qt, we may replace C by C G L^(0, T; L^(12)) in (51), where 
C is given by 

C = -X{z=o}max|o, W^(e(u),z) -k /'( 2 )|. 

We check that C satisfies (13a), i.e., C is a desired subgradient. 


j ^VZr(t) • VCr(t) + {W^zisiUrit)), Zr{t)) + f'{Zr{t)) -k 5t£r(t))Cr(i)) 


dx 
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To (27): Let ti and t 2 with 0 <ti < t 2 <T he arbitrary. Integrating (47) from Lemma 3.15 
over the time interval [ti,t 2 ] yields 



^7>(t) + + 2 ?t-(0, t) + Tr(0, t)^ dt < 



(.T(o)+/c(o) + >v;,i(o,t)) 


dt 


By the convergence properties in Lemma 3.12 and by lower semi-continuity arguments, 
we obtain 


rt2 

liminf / 


r \0 


'ti 


/*^2 pt'Qi 

(t) dt + liriMnf J ICr{t)dt> J ^F{t) + dt. (52) 


The limit passage in the dissipation term f Pr(0, t) can be performed by Fatou’s lemma, 
the estimate tr > t, the weak convergence dtZr dtz in L^(0,T;L^(ll)) (see Lemma 
3.12) and by a lower semi-continuity argument: 


rt2 

lim inf / V. 


(0, t) dt > liminf / / / |9t2'T('S)P dxds dt 

Jti Jo Jn 

liminf [ [ [ 1942^(5)1^ dxdsdt 

Jti Jo Jn 

f (^liminf f f |(9tZr(s)p dxds^ dt 

Jti V Jo Jn J 

f f f |5tz(s)p dxdsdt. 

Jti Jo Jn 


> 


> 


(53) 


(54) 


Moreover, Lemma 3.12 and Lebegue’s convergence theorem lead to 

hm r>VI,,(0,t) = 0. 

^\o Jti 

To treat the error term f Sr(0,t}, we define Tr(0,t) =: S^(0,t) +S^(0,t) with 

£^(0,t) := f f - ^ ^ ^ -^-^Ce(tt“) : e(n7) dx ds-|- f f W^z{£{ur), Zr) dtZr dx ds 

Jo Jn^ 'T Jo Jn 


T 

- r / L£lLT(iLdi.ds+ 

Jo Jn T 


f ^ [ f'i^r) 

Jo Jn 


dtZr dx ds. 


By the differentiability of /i, it holds 


h{z.^ ) = h{Zr) + h'{Zr){z.^ - Zr) + r{z.^ - Zr), 


r{r]) 


0 as T/ ^ 0. 


We then get 

'' 1 h{z~) - h{zr) 


r f -■ 

Jo Jn 2 


r 


/7 

Jo Jl; 


Ce(tt.,-) : £{u .^) dx ds 

1 


{z^ (s)7^Z.r(s)} 


h'{Zr) 


Z.r- — Zr r[Zr — Zr) Zr — Z- 


+ 


Zr — Zr 


C£(u .^) : e{u .^) dx ds 
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(55) 




Zr)— - -C£{u^ ) : £{u^ ) dx ds 

T 


=.Ti 

1 r{z~ - Zr) Z~ - Zt 
/o J {z^ {s)^Zt{s)} 2 Zr Zr 


+ 


/7 

Jo J{z 


C£{u^ ) : £{u^ ) dx ds 


= -.T2 


By using the convergence properties in Lemma 3.12, we find 


Ti 


f f -h'{z)C£{u) : £{u)dtz dx ds 
Jo Jn 2 

= f f W^z{^{u), z)dtz dx ds 

Jo Jn 


as T \ 0. 

The Lipschitz continuity of h on the interval [0,1] implies the boundedness of 


riZr - Zr) 


Zr — Zr 


< 


L°°{{Zr ^Zt}) 


h{Zr ) - h{Zr) 


Zq- Zq- 


+ 


L°°{{Zr ^Zt}) 


h'{Zr)-^ 


Zr — Zr 


L°°({2r ¥^^r}) 


< c. 


Taking also —)• 0 a.e. in 0 x (0, T) as r \ 0 into account, we conclude by 

Lebesgue’s generalized convergence theorem 


r{Zr - Zr) 


Zr — Zr 


Ll{{Zr ^Zr}) 


0 for every q G [l,oo). 


(56) 


Therefore, we find by Holder’s inequality 

\r{Zr{s) - Zris)) 


T2< 


-I 

2 Jo 


Zr (s) - Zris) 


dtZr{s)\\L2(^n)\\£{Ur (•?)) ||L4(f1) 


< c 


r{Zr - Zr) 


Zr — Zr 


L4({2^ ^2r}) 


L'^{{Zr (s)^2r(s)}) 

*t^T||L2(o^r;L2(Q))lkr llL°°(0,T;H2(f2;R")) 


7r"^ IL-»({ 2 -^ 2 r}) ^0^1^ (^6) aud 11<9tZr 11 1 , 2 (o,t;L 2 (Q)) as well as ||u^ ||L°c(o,T;H2(n;R-)) 
are bounded by Lemma 3.10, we obtain T 2 —)• 0 as r \ 0. The convergence properties 
in Lemma 3.12 also yield 



W^zi^iUr), Zr) dtZr dx ds —)• 


[ [ W^z{e{u),z)dtz 

Jo Jn 


dx ds. 


In particular, we have used £{ur) £{u) in L^(0, T; L'^(n; M"^"')) due to the H"^- 

regularization. 

Together with (55) and the identified limits for the terms Ti and T 2 , prove £r{0, t) —)• 0 
as r \ 0 for all t G [0,T]. Taking also the uniform boundedness |T.,1(0, t)| < C with 
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respect to r > 0 and t E [0, T] into account (which follows from the a priori estimates 
in Lemma 3.10), we find 


j•t2 

lim / t) dt = 0 

Jti 

by Lebesgue’s convergence theorem. Note that the limit 

rt2 


rt2 

/ £Uo,t) 

J t1 


dt = 0 


(57) 


(58) 


lim 

can be shown with the same arguments. 

Integrating (47) over the time interval [ti,t 2 ] and applying liminf^-^o both sides, we 
obtain 

rt2 rt2 rt2 rt2 

lim inf / J>(t) dt + liminf / /CT-(t) dt + liminf / Pt-( 0, t) dt + liminf / £’t-( 0, t)dt 

Jti Jti Jti Jti 


< 


rt2 




+ 




dt + lim inf / A 

t\0 Jti 


The convergence properties (52), (53), (54), (57) and (58) show 


r^2 


^2 


(.F(t) + /C(t)+P,(0,t))dt < (.T(0) + /C(0)+>Vext(0,t))dt. 

Since ti and t 2 with 0 <ti <t 2 <T are arbitrary, we obtain (27). 

Hence, we have established existence of weak solutions in the sense of Definition 3.7. 


□ 


3.3 Existence of weak solutions for the limit system 

In this section, we are going to prove Theorem 3.6 by performing a limit analysis 5 0 

for solution {us,zs,^s) from Proposition 3.8. To this end, we approximate the data {vP,i,b) 
given in Theorem 3.6 by smooth functions (e.g. via convolution) G is G 

C°’i(0,r;L2(D;M^)), bs G ^^^(O,T; 772 (^ 7 -such that 

—)• strongly in 77^ (D;M"'), (59a) 

is ^ i strongly in L‘^{0,T; (59b) 

bs^b strongly in 77^(0,r;77^(H;M”)n 77^(0,r;L2(D;M")) (59c) 

as (5 0. By possibly reparametrizing we obtain the following a priori estimate 

uniformly in 6: 


V~6\\ul\\Hr < a (60) 

The cornerstone of the passage 5 0 in the weak formulation are the following a priori 

estimates for us and zs uniformly in 6 which are obtained by means of the regularized total 
energy-dissipation inequality (27). 
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Lemma 3.16 (A priori estimates) The following a priori estimates hold uniformly in 5: 


V^lk<5||L°°(0,T;/f2(n.iRn)) < C, (61a) 

°°(0,T;/fl(C;R"))nVKl.°°(0,T;L2(0;Rn))p|j^2(o,T;(//2^(0;R"))*) — C*, (61b) 

lk5||L°°(0,T;VKl'P(0))nHi(0,T;L2(n)) < C. (61c) 


Proof. By the a priori estimates in Lemma 3.10 (i) which are independent of r and 6, by 
lower semi-continuity of the norm and by (60), we find (Ola) and (61b). 

To gain the a priori estimate (61c), we use the total energy-dissipation inequality (27). It 
remains to estimate the following terms occurring in the Wf^-j-term: 


n W^e(£^{us), zs) : s{dtb)dxds < ||ti5||L2(o,T;/fi(f^;i*’^)) 

! 

n dtU 5 ■ dttbdxds < \\us\\H^(^o^T;L'^{n;K'^))\\b\\H2{0,T-,L^{n-,R^))^ 

! 

/ dtusit) ■ dtb{t)dx < ||u5||vi/i’°°(o,r;L2(r2;R'i))||^||vi/i-°°(o,T;i,2(n;Rn)), 

Jn 

n h ■ dtius — b) dxds < \\h\\L^{o,T-,L^{n-,R’^))\Ws — 

! 

6 / {Aus{t),dtb{t))jj 2 dt < (5||t(5(t)||^oo(o,T;//2(Q.Rn))||9t6||^oo(o_T;H2(0;R'*))- 

Jo 

Taking the estimates (61a) and (61b) into account, we see that t) is uniformly bound¬ 
ed in t and 5. Therefore, by (27), J^s{t) is also uniformly bounded which proves (61c). □ 


Lemma 3.17 There exist functions 

u e L°°(0, T; (0; M")) n IT^’°°(0, T; L^{TL- M")) n 77^(0, T; (77^^ (17; M”))*), 
z G L°°(0, T- IT^’P(I7)) n 77^(0, T; 7.2(17)) 

satisfying (9), (12b), (13b) and a subsequence {6k}keN with 5k\0 as k co such that 

us^ weakly-star in L°°{d,T]H^{Tt]W^))r\W^’°°{d,T]L^{n]MT)) 

and weakly in 77^(0, T; (77^^(17; M"))*), (62a) 

zsf. ^ z weakly-star in L°°{0,T;W^’^{n)) and weakly in H^{d,T]L^{Tt)), (62b) 

^4 ^ z strongly in 7.^(0, T; IT^’P(f7)), (62c) 

zs^. z strongly in C(fl x [0,T]) (62d) 

as k y' oo. 

Proof. Properties (62a), (62b) and (62d) for some functions u and z follow with standard 
compactness and Aubin-Lions type results (cf. [Sim86]) by keeping p G (n, oo) in mind. 

We also obtain by Aubin-Lions 

zsf. z strongly in L^(0,T;7.°°(f7)) for every q G [l,oo) (63) 
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and thus (by possibly resorting to a further subsequence) 

Z 5 f,{t) ^ z{t) weakly in (64) 

as k Z' oo and for a.e. t S (0,T). 

The strong convergence (62c) can be shown as in Lemma 3.14: We apply Lemma 2.1. Due 
to (64), one finds an approximation sequence {C 4 } ^ L°°( 0 , T; W:J^’^(D)) satisfying Csk C 
in L^(0, T; 1 T]l’^(D)) and 0 < C 4 < -24 a.e. in D x (0, T). By using the same p-monotonicity 
estimate as in (45), testing (21) with ~ -2^4) using the a priori estimate (61b) in Lemma 
3.16 and the convergence (63), we obtain (62c). □ 


As usual, we omit the subscript k. 

Remark 3.18 Note that by (28) and (61b), we only obtain an (0,T] (^}))-bound for ■ 
A major challenge in the passage 6 is to establish a desired subgradient G L^(D x (0, T)) 
for the limit system. 


Proof of Theorem 3.6 

We are going to prove that the functions u and z from Lemma 3.17 also satisfy (10), (13a), 
(21) and (23). 

To (10); Integrating (26) over time from 0 to T and using the definition for the elastic energy 
density in (4), we find 


[ {dttus{t),Cit))H^ df+ f f h{z 5 )Ce{us) : e{C) dx-\-5 f [ (V(Vn 5 ), V(VC)) dxdf 
Jo D Jo Jn Jo Jn 


[ [ 4-Cdx 
Jo Jn 


for all C £ By exploiting the convergences (59b), (62a), (62d) 

and 


<5 



{V{Vus),V{VC))dxdt 


< <5|k5||L2(0,T;/i'2(f7.]Rnp ||C||L2(o,r;i^2(0;M<i)) —> 0 


due to (61a), we conclude (10) for all f G and a.e. t G (0,T). 

By using the density of the set L7p^(D;M”) in L7p^(D;M") (here we need the as¬ 
sumption that the boundary parts To and Tn have finitely many path-connected 
components, see [Berll]), we identify dttu{t) G (L7p^(D; M"))*. Consequently, the 
equation (10) is true for all C, G (D;M") and a.e. t G (0, T). In particular, 
a«nGL°°(0,r;(Ri^(D;M-))*). 

To (13a) and (21): We choose the following cluster points 


X5 X{zs>0} X 

Vs •= X{zs=0}n{W,z{£{us),zs)+f'{^s)<0} V 

Fs ■■= X{zs>o}^ius) F 

Gs ■= X{z5=o}n{M/.,(£K),2i)+/'(25 )<o}^('«<5) G 


weakly-star in L°°(D'r), 
weakly-star in L°°(D 7 ’), 
weakly in L^(D 7 ’; M”'^"'), 
weakly in 


(65) 

( 66 ) 

(67) 

( 68 ) 
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as (5 \ 0 for a subsequence. By (62d) and (62a), we obtain for a.e. x £ {z > 0} 


Xix) = 1, r]{x) = 0, F{x) = e{u){x), G{x) = 0 (69) 

by the following comparison argument: 

Let C £ with supp(C) F {z > 0}. Then, by (62d), we find supp(C) ^ 

{zs > 0} for all sufficiently small <5 > 0. On the one hand, (67) implies 

f f : C dx dt —)• f [ F : ( dx dt. 

Jo Jn Jo Jn 

On the other hand, by supp((^) C {zs > 0} and (62a) 

f f Fs:Cdxdt= f f e{us):(dxdt^ ( f e{u):Cdxdt. 

Jo Jn Jo Jn Jo Jn 

Thus Jq f^^iu) ■ Cdxdt = JqF : Cdxdt and, consequently, F = £(u) a.e. in 
{z > 0}. The other identities in (69) follow analogously. 

Now, let C G L°°(0,T;W^’^(Q)}. Taking (28) into account, inequality (21) becomes 
after integration over time 

0< / [ {\V Z5\P~‘^V zs ■ V C + dtzsC) dx dt + [ {W^z{£{us),zs) + f'{zs)) (dxdt 

Jo Jn J{zs>o} 


+ 


/ {W^:,{e{us),zs) + f'{zs)) Cdxdt (70) 

J {zs=0}r\{W,z{e{us),zs)+f'{zs)<0} 

for all C £ L°°(0, T; lyl’^(n)). Applying limsnp^^o on both sides and mnltiplying by 
— 1 yield 

0 > lim /" [ {\Vzs\^~‘^Vzs ■ V{-() + dtZs{-C)) dx dt 

s\oJo Jn 

Tliminf [ [ h'{zs)CFs : Fs{-C) dx dt + Vun [ [ xs f{zs){-C)dxdt 


5\o Jo Jn 

+ lim inf / / 

^\o Jq Jq 


s\oJo Jn 

rT 


h'{zs)CGs : Gs{-C)dxdt + l[m f f rjs f'{zs)i-C) dx dt. 

^\0Jq Jn 


Weakly lower semi-continuity arguments, the convergence property (62d) and the iden¬ 
tifications listed in (69) give 


0 > 


>0 Jn 


{\Vz\P-Wz • V(-C) + dtz{-C)) dxdt 


+ 


>{z>0} 


(W^(e(n), z) -I- f'{z)) (-C) dx dt 


+ [ (h'{z){CF ■.F + CG:G) + {x + v)fiz)) {-() dx dt. 

J{z=0} ^ ^ 


This inequality may also be rewritten in the following form: 


+ / (h'{z){CF :F + CG:G) + ix + v)f'iz) - W,Mn),z) - f{z)) C dx dt. 

J{z=0} ^ ^ 

Therefore, 


j (l Vz|^“Vz • VC + {W,,ieiu),z) + fiz) + dtz + C) C) dx dt 


with 


C := X{^=o}min|o,/i'(2;)(CF : F + CG : G) + {x + V - l)/'(2) - VT,^(e(u), 2;)|. 

We obtain (13a) and (21). 

To (23): In order to establish the energy inequality (23), we can proceed as in the proof of 
Proposition 3.8: 

In fact, integrating the regularized energy-dissipation inequality (27) over the time 
interval [ti,t 2 ] with 0 <ti <t 2 <T, using the estimate ^{Aus{t), us{t))jj 2 > 0 on the 
left hand side and passing to the limit by using lower semi-continuity arguments, the 
convergences in Lemma 3.17 and (59a)-(59c) as well as Fatou’s lemma yields 

r {F{t) + }C{t) + P(0, t)) dt < r {F{d) + /C(0) + We.t{0, t)) dt. 

Jtl Jtl 

Since 0 <ti <t 2 <T were arbitrary, (23) holds for a.e. t £ (0, T). □ 
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